Interaction induced collapse of a section of the Fermi sea in in the zig-zag Hubbard 

ladder 
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Using the next-nearest neighbor (zig-zag) Hubbard chain 
as an one dimemensional model, we investigate the influ- 
ence of interactions on the position of the Fermi wavevec- 
tors with the density-matrix renormalization-group technique 
(DMRG). For suitable choices of the hopping parameters we 
observe that electron-electron correlations induce very differ- 
ent renormalizations for the two different Fermi wavevectors, 
which ultimately lead to a complete destruction of one section 
of the Fermi sea in a quantum critical point. 



Introduction - Many aspects of the low-energy 
physics of an electronic system are influenced by the 
shape of its Fermi surface and the occupation of nearby 
states. Associated observables may vary strongly with 
the temperature and the strength of the interaction of 
the electrons. Changes in the Fermi-surface geometry in- 
duce magnetic and other instabilities. Electron-electron 
interactions induce a momentum dependent softening of 
the Fermi surface, which is responsible for a wide vari- 
ety of phenomena, such as the loss of magnetic order or 
high-temperature superconductivity. In systems where 
the Fermi surface is not simply connected, interactions 
may lead to a partial or total collapse of parts of the 
Fermi surface at a quantum critical point and an asso- 
ciated drastic change in the physical properties of the 
system. It is therefore important to study the renormal- 
ization of the individual Fermi-surface sections under the 
influence of electron-electron correlations in competition 
with frustrating interactions. 

Recent investigations suggested a spontaneous, inter- 
action induced deformation of the Fermi surface of the 
of the 2D t-J model Q] as well as of the (extended) 
Hubbard model |^-|| . Because the effect of strong elec- 
tronic interactions is difficult to study in two- or three- 
dimensional systems, these studies are confined either to 
the weak coupling limit or involve extensive numerical 
studies of effective lattice models. Many aspects regard- 
ing the Fermi surface renormalization for strong interac- 
tions remain presently unclear. 

The development of the density matrix renormaliza- 
tion group now offers a reliable — though numerically 
involved — approach to rigorously study the effect of 
strong interaction in one dimension. In this letter we re- 



port a study of the renormalization of the Fermi wavevec- 
tor, or Fermi point, in the Hubbard chain with next- 
nearest-neighbor interactions as a one-dimensional model 
for Fermi-surface renormalization. We confirm results 
of a recent renormalization group (RG) analysis for the 
Fermi point renormalization in weak coupling. Fermi 
points close to a saddle-point of the momentum distri- 
bution function n(k) are predicted to shift towards the 
van-Hove singularity. A similar renormalization has been 
suggested for two-dimensional systems jlJQ]. For strong 
interactions, the pocket of the Fermi-sea near the saddle- 
point may be emptied completely, a hypothesis Q which 
can not be verified in weak-coupling. In order to study 
this prediction rigorously, we have performed a DMRG- 
study of the t\ — ta Hubbard model for rings with pe- 
riodic boundary conditions, evaluating the momentum 
distribution function n{k) directly for rings with up to 
80 sites. For larger interaction strength we find a novel 
interaction-induced quantum critical point that is asso- 
ciated with the collapse of a Fermi sea pocket. This phe- 
nomenon is analogous to the opening of a pseudogap, 
found experimentally in the high-Tc-cuprates @, near 
the saddlepoint of the electronic dispersion. 

The Hamiltonian of the zig-zag Hubbard ladder (see 
Fig. (fil), top) is given as: 
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where the a (c n CT ) are Fermion-creation (destruction) 
operators on site n and spin a =\,i. For appropriate 
choices of the parameters this model describes either the 
low-energy properties of Hubbard-ladders |(| or half-filled 
edge-sharing double-chain materials like SrCu02 § or 
LiVaOs 0] for which the next-nearest neighbor hopping 
ty, is expected to be substantial. 

Generically, the variation of the two hopping param- 
eters (£1,2) allows changes in the position of the Fermi 
point in the non-interacting limit. When the hopping t-i 
is large enough, the Fermi-sea will split into two separate 
parts, see Fig. ([!]). According to a recent RG analysis for 
the case of two separate Fermi-seas Q the Fermi points 
kp (i = 1, 2) are renormalized as: 
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In weak coupling, the Fermi point-shift Ak depends only 
on the respective Fermi-velocities Vi and on the initial 
momentum cut-off Ao. Eq. (||) is rigorous in the weak 
coupling limit. 

DMRG - The evaluation of n(k) as the Fourier trans- 
form of the correlation function: 



n a (k) 
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is numerically difficult and costly in the framework of the 
DMRG Jl^,^l|, in particular for periodic boundary con- 
ditions |p.2| . For large interaction strength, the suppres- 
sion of double occupancy reduces the size of the relevant 
Hilbcrt space and improves the convergence properties of 
the DMRG procedure |l3| . We have therefore decided to 
study the case of quarter-filling, where the renormaliza- 
tion of the Fermi point occurs at relatively large values 
of the interaction U. We chose t 2 /ti = —3 as a good set 
of parameters to observe the Fermi point renormaliza- 
tion effects at quarter-filling, since one of the Fermi-seas 
is small for this case (see Fig. ([[])) and since the weak- 
coupling RG predicts gapless Luttinger-liquid behavior 
for these parameters B . 
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FIG. 1. Top: Illustration of the t\-ti zig-zag lad- 
der, the dashed/full lines correspond to t\/t2 bonds. Bot- 
tom: Illustration of the U — dispersion relation 
E(k) = 2ti cos(fc) +2t2 cos(2fc) for the case ti = 1, ti — —3. 
Also show (horizontal dashed line) is the Fermi-energy 
Ef = — 3.6736ti for the quarter-filled case. 

We checked for both the convergence of the variational 
energy and of the correlation functions by increasing the 
number of DMRG states m and found reasonable choices 
for the desired accuracy. The weight of the discarded 
fraction of the density matrix varies between 1 x 10~ 4 to 
2 x 10 -5 depending on the choice of system size and the 
number of states. Fig. (0) demonstrates the convergence 



of n(k) both as a function of lattice size and of the num- 
ber of states m in the DMRG procedure. Additionally, 
we performed calculations for open chains and verified 
that the behaviors of the real-space correlation function 
is consistent with those for periodic systems. 
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FIG. 2. Convergence of the correlation function n(k) com- 
pared for different chain lengths L and DMRG states m (in- 
set shows L — 80) for a quarter-filled system with t\ — 1, 
t2 — —3. The lines are guides to the eye. 

Results - Depending on the relative strength of the 
hopping parameters £2/^1 and interaction strength U/t\ 
the ground state phase diagram of the t\ — <2-modcl 
exhibits various instabilities towards the formation of 
states with spin- and/or charge-gap JT^JT^] . We have 
therefore investigated the possibility of instabilities to- 
wards ferromagnetism and phase-separation by evalu- 
ating AE F = E(N 1 + 1,N X - 1) - E(Ni,N{), where 
E(N<t, Ni) is the ground-state energy of the quarter-filled 
system, and AE ph = E(N 1 + 1, JVj. + 1) + ^(A^t ~ 1 , N l~ 
1) — 2E(N^, Ni). For the parameters considered here no 
tendency towards an instability was found for any inter- 
action strength, in accordance with Ref. fl5|| . 

In Fig. ([}]) we present results for the momentum dis- 
tribution function for U jt\ = 1,30,100. For ti = — 3£i, 
quarter filling and {7 = the two Fermi wave-vectors 
are kp — 0.055-7T and kp — 0.8057T. The parame- 
ter a = (vp + v^)/(2vjp) — 5.53 entering the weak- 
coupling RG-equations M indicates a gapless Luttinger- 
liquid phase at small couplings. 

For small U we clearly observe the existence of 



two Fermi-points kp^ and kp (see Fig. 3), which 
are substantially renormalized with growing interaction 
strength. 

It is evident from the data presented in Fig. (|^), 
that the magnitude of the momentum distribution func- 
tion inside the smaller Fermi-sea is reduced substantially 
stronger by the interaction than the one inside the large 
Fermi-point. This is a consequence of the very different 
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values for the respective Fermi-velocities, v\ — 0.658 1\ 
and vi = 6.620 ti, which differ by one order of magnitude. 
The energies of particle-hole excitations are ~ vpAk, 
many more particle-hole excitations are therefore created 
for the Fermi-sea pocket, which has the smaller vp = V\, 
resulting in a substantial reduction in the occupation 
numbers n(k) for k < kp . 

In order to quantify the results of Fig. 3 we analyze the 
Fermi point position as the point where the slope of the 
momentum-distribution function is maximal. We find 
that our finite-lattice data for the momentum distribu- 
tion function n(k) is approximated well by two washed- 
out step-functions: 
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The quality of this fit is illustrated by the lines in Fig. 
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FIG. 3. The DMRG-results (symbols) for the momentum 
distribution function for a quarter-filled 80-site system with 
periodic boundary conditions. The parameters are ti = 1, 
ta = —3, and U — 1 and U = 30 respectively. The lines 
through the DMRG-data are fits by Eq. (Q). The dashed verti- 
cal lines indicate the (7 = positions of the two Fermi-points. 

Fig. (Q) shows the dependence of kp^ 2 ^ as well as the 
total volume of the Fermi SGet clS c\i function of U (in 
units of 7r). The data clearly indicate the presence of 
a quantum-critical-point with U c ~ 50ii, for which the 
smaller Fermi sea collapses. The fit-parameter kp enter- 
ing (|[) denotes a locus of maximal slope in n(k), which 
corresponds to the Fermi-point position for U < U c . As 
the calculations were done for fixed particle density, we 
expect the Fermi-sea volume to be independent of the 
interaction-strength U (Luttinger's theorem |l6|). From 
Fig. (||) we see that this expectation is approximatively 



fulfilled if we take (tt — kffi) + kp' for the volume of the 
Fermi-sea for U < U c (open diamonds in Fig. (0)) and 

(2) . . 

and tt — kp for U > U c . Some numerical uncertainty 
is observed, as usual in numerical simulations, near the 
critical point. The error bar in the figure denotes the lat- 
tice induced finite-size error on the Fermi point position. 
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Within this error all data is compatible with the exact 
result in the noninteracting limit. 

As the data for U = 100 in Fig. 3 indicates, we have 
also maxima in the slope of n(k) corresponding to the 
remainder of the small Fermi-sea for U > U c . To distin- 
guish between incoherent excitations and a sharp Fermi 
point we have plotted the hypothetical total Fermi sea 
volume (open circles in Fig. 4) under the assumption that 
these inflection points also correspond to Fermi points. 
As the figure indicates this would lead to a massive vio- 
lation of the Luttinger theorem for all U > U c . We thus 
conclude that n(k) for small k and U > U c can be con- 
sistently associated with contributions from incoherent 
excitations to the momentum distribution function. 
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FIG. 4. The DMRG-results for the volume of the small 
(fc^ 1 '. filled squares, the line is a guide to the eye) and large 
(tt — kp" 1 , filled diamonds) Fermi-sea and the total (open dia- 
monds) volume of the Fermi-sea (FSV) . The latter is identical 
with the volume of the remaining Fermi sea for U > 50. The 
dashed line indicates the analytical value of the total FSV 
at U — (quarter filling). The open circles indicate the total 
FSV under the assumption that a a small Fermi point persists 
for U > 50. 



Conclusions - In a one-dimensional model we have 
analyzed the stability of Fermi-sea pockets with respect 
to electron-electron interactions and found that suffi- 
ciently strong repulsive interactions lead to a novel quan- 
tum critical point at which the one Fermi-sea pocket is 
destroyed altogether. We studied this quantum-phase 
transition within the t\ — ti zig-zag ladder as a proto- 
typical model. We estimate the value of U c /t\ » 50 for 
t-z = 3<i in a quarter-filled system. We note that the 
nature and existence of this transition does not depend 
on the specific choice of 4i,i2 an d the filling, but that 
a manifold of transitions exist for appropriate choices of 
these parameters. The conditions of the existence of the 
transition are the existence of two different pockets of 
the Fermi sea with a significantly higher Fermi velocity 
at the Fermi point of the smaller pocket. In our study we 
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have chosen the values of the parameters such that the 
transition becomes amenable to quantitative treatment 
with the DMRG, which converges progressively worse 
near half-filled systems. This lead to a critical U c that 
seems to be too large to be relevant for experimental 
realizations of (Q), but we would like to point out that 
at least two effects are likely to dramatically reduce the 
critical U c : 

(i) At half filling the £2 — ti model undergoes a Mott- 
Hubbard transition for ti < ti/2. For ti — > ti/2 the 
critical U c — > 0. It has been estimated || that the 
RG-estimate for the critical U c obtained from Eq. (|^) [12 
is in agreement with numerical results 0Jl^] . We there- 
fore expect the critical U c for the destruction of the 
small Fermi-surface pocket to reduce drastically near 
half-filling. Here we did not investigate this region due 
to numerical difficulties. 

(ii) In realistic systems the Coulomb-interaction will 
be longer-ranged than the onsite-form assumed in the 
present study. It is known, that longer-range contribu- 
tions to the interaction increase substantially the effect of 
the interaction on the Luttinger-liquid parameters in the 
one-dimensional electron gas |l9|,^0|]. We expect a simi- 
lar enhanced influence on the renormalization position of 
the Fermi points. 

We believe our rigorous findings for the one- 
dimensional case to be relevant for two-dimensions. We 
note that the Fermi-surface instability observed in the 
present study is qualitatively different from the Pomer- 
anchuk instability observed in weak-coupling |2lJ] for the 
2D Hubbard model, since the Pomeranchuk instability 
corresponds to a spontaneous breaking of the tetrago- 
nal (C4) symmetry due to enhanced scattering from one 
saddle-point to another. The Fermi-surface instability 
observed in the present study is driven, on the other 
hand, by scattering between 'normal' and flat parts of 
the Fermi-surface. 
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